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Abstract 

An effective equation of state which generalizes the Lattimer-Swesty equation for 
nuclear matter is presented for matter at supernuclear densities including strange 
baryons. It contains an adjustable baryon potential energy density, based on models of 
local potentials for the baryon-baryon interactions. The features of the equation rely on 
the properties of nuclei for the nucleon-nucleon interactions, and mainly on experimen- 
tal data from hypernuclei for the hyperon-nucleon and hyperon-hyperon interactions. 
The equation is used to calculate equilibrium compositions and thermodynamic prop- 
erties of high density matter with strangeness in two astrophysical contexts: neutron 
star matter (transparent to neutrinos) and proto-neutron star matter (opaque to neu- 
trinos). The effective equation of state reproduces typical properties of high density 
matter found in theoretical microscopic models. Of these, the main result is that hy- 
perons appear in both types of matter at about twice the nuclear saturation density, 
and that their appearance significantly softens the equation of state. The range of 
maximal masses of neutron stars found in a comprehensive parameter survey is 1.4-1.7 
Mq. Another typical result is that the maximal mass of a proto-neutron star with 
strange baryons is higher than that of an evolved neutron star (opposite to the case of 
nuclear matter), setting the stage for a "delayed collapse" scenario. 
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1 Introduction 



The properties of matter at supernuclear densities are a key subject of interest in the study 
of core-collapse supernovae and neutron stars. Theories of such matter typically predict that 
at the nuclear saturation density, ps = 0.15 — 0.16 fm~^, it is composed of free nucleons and 
leptons, but at higher densities several more species of particles may appear. Chiefly among 
these are strange baryons, namely, the A, S and S hyperons, while A baryons along with 
pion and kaon condensations and a deconfined-quark-phase have also been considered. The 
collapsed core in type II supernovae and neutron stars serve, in this sense, as cosmological 
laboratories for hadronic physics. 

The appearance and effects of such "novel" particles in dense matter have been included in 
several studies in the context of neutron star structure (see for example, Refs. [1-5]). Some 
work has also been done regarding the collapsing core in type II supernovae , and recently 
attention has also been given to the possible impact of the appearance of strange particles 
in proto-neutron star evolution [7-10]. 

The physical state and composition of matter at high densities depends mainly on the na- 
ture of the strong interactions. However, the theoretical details of these interactions are not 
completely understood, and uncertainties also arise from the incompleteness of the available 
experimental data for nuclear matter. Such uncertainties are compounded when strange 
hadrons are taken into account, as relevant experimental data is rather scarce. Theoreti- 
cal studies of nuclear and hadronic matter are usually performed with sophisticated models 
such as relativistic mean field (RMF) and Brueckner-Hartree-Fock, which involve effective 
Lagrangians and/or specific assumptions about baryon-baryon interactions. Some applica- 
tions of these models are technically complex and often prove expensive in terms of CPU 
time when performing calculations at various combinations of densities and temperatures. 
Thus, there is also an interest in using an effective equation of state (EoS), which is ana- 
lytic and depends on a small number of parameters. Such an equation will be convenient 
as a supplementary tool, mainly for direct use in rapid hydrodynamical simulations, and for 
conducting extensive parameter surveys that are necessary given the large uncertainties in 
the theory and experimental data regarding strong interactions. 

Two commonly used effective equations of state for hot (finite temperature) nuclear matter 
(nucleons and leptons) are those of Lattimer and Swesty |ll|] and Baron, Cooperstein and 
Kahana [|T^]. In the present study we wish to construct a generalized effective equation of 
state for high density hadronic matter, by extending the Lattimer-Swesty (hereafter LS) 
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equation so as to include additional species of particles. We currently limit ourselves to 
strange baryons, which we view as the prime candidates for appearing in high density matter. 

In Section 2 we discuss our approach for constructing the effective EoS, basing it on the for- 
mulation of local (effective) density dependent potentials for the various cases of the baryon- 
baryon interactions (lepton interactions are assumed to be negligible, and thus leptons are 
treated as non- interacting particles). Section 3 outlines the equilibrium conditions, which 
determine the composition and physical properties of high density matter for neutrino-free 
matter (neutron star) and neutrino trapped matter (core-collapse and proto-neutron-star) 
in beta equilibrium. 

Our choice of parameters for these equations, derived mainly on the basis of hypernuclear 
data, is presented in Section 4; we note, that these parameters allow accurate modeling of 
field-theoretical interactions. We also specifically treat the consequences of recent theoreti- 
cal and experimental work on S~-atoms [|l^, [l^, which seem to suggest a strong isoscalar 
repulsion in the interaction of S hyperons with a bulk of nucleons, unlike the commonly 
assumed interaction. 

Results for the equilibrium compositions of high density matter and the corresponding EoS 
are presented in Section 5. These EoS's are examined with regard to the static neutron star 
masses they predict in Section 6. In Section 7 we offer conclusions and further discussion. 



2 The Effective Equation of State 

The construction of an effective EoS for supernuclear density matter centers on the descrip- 
tion of the strong interaction energy. In this approach one models these interactions in some 
analytic form, involving a small number of parameters. 

The internal energy density is approximated by dividing it into a baryon contribution, which 
includes kinetic, potential and mass terms, and a lepton contribution which has only kinetic 
and mass terms. This approximation, introduced in the LS EoS, yields the following structure 
of the internal energy density for supernuclear dense matter: 

£tot = £kin{p, T, {Xi}) + Epotip, {Xi}) + Emassip, {Xi}) + Ekin{p, T, {k}) + Emass{p, {k}) , (1) 

where p is the total baryon density, T is the matter temperature (whose effect is handled 
in the kinetic energy terms), {xj} and {/j} are respectively the baryon and lepton fractions, 
Xi = Px J P] h = Pk/ P- For finite temperatures, a photon contribution must be added as well. 

The separation of the lepton contribution is made possible by treating leptons as non- 
interacting (relativistic) particles, and, as discussed in Ref. [l^, electron screening effects 
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may be ignored because the electron-screening length is larger than the baryon separation 
distance. 



The heart of the effective EoS is, of course, casting the potential energy density term. Spot, 
into an effective form. This is achieved by assuming local density dependent potentials for all 
strong interactions. Such a potential model must reproduce the basic features of the strong 
interactions, i.e. long-range attraction and short range repulsion, and - in some cases - charge 
dependence, compatible with isospin invariance. The most economic way to construct such 
a potential for the case of a single baryon of species y in bulk matter of baryon species x of 
number density is: 

OjxyPx ~\~ ^xy^x^yPx ~\~ ^xyPlc ^ ■ (2) 

The first term yields attraction [a^y negative), the third term yields repulsion [cxy positive, 
supposedly introducing multibody interactions), and the second (symmetry) term introduces 
the charge dependence through a charge (isospin) t. The common choice of parameters 



in effective approaches |I2| favors a linear density dependence of the attraction and 
symmetry terms, while 'jxy must be greater than unity, so that repulsion will dominate 
at high densities (short ranges). This formulation limits the potential energy model to at 
most four parameters, whose values are chosen to reproduce experimental data and accepted 
theoretical results. Eq. (|^) can, of course, also be used for the case of a single species (i.e., 
y = x), describing the local potential felt by each baryon in the bulk of baryons of the same 
species. 

We note that a local effective potential for the nucleon-nucleon interaction provides the basis 
for the well known Skyrme model (see Ref. [jl6l for a recent review). Similar modeling for 



hyperons in nuclear matter was first performed by Millener, Dover and Gal ([|T7|, hereafter 
MDG), on the basis of experimental data of various A-hypernuclei, deriving a local potential 
for a single A in nuclear matter. 

The local potential for a single baryon in a bulk of other baryons may be extrapolated into 
the potential energy density of bulk matter with a total density p that includes both types 
of baryons. This is done by folding Vy{px) with the partial density Py, and vice versa (K(Py) 
with px), combined with weight factors (avoiding double counting of each interaction). If 
we choose the weight factors to be x/(x + y) and y/{x + y) respectively, the contribution of 
the interactions between baryons of species x and baryons of species y to the internal energy 
density, e{x ^ y,p) will be: 

£pot{x ^y,p) = axyxyp^ + bxytxtyXyp^ + Cxy i^j^r-^^'^'''' V + '^^~^y^''''"^ pT-iz+i , (3) 
where x and y are defined as x = px/ p, y = Py/ p. 
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The above weight factors do not appear in the first two terms as they just add to unity, while 
in the repulsion term they retain discrimination between baryons y in matter of baryons x 
and baryons x in matter of baryons y. Although this form is not unique, the weight factors 
x/{x + y) and y/{x + y) yield the expected results in the case x = y and for the limits 
x»y and vice versa. Inherent in Eq. (^) is the neglect of possible screening effects of other 
baryons, so that the interaction between any two baryons is dependent only on the distance 
between them. This interaction may be taken as an effective average interaction between 
any two baryon species. 

For the specific case of baryons of a single species, setting x = y m. Eq. (|^) requires a factor 
of \ to avoid double counting: 

epot{x ^ x, p) = ^ a^xX^p^ + h^^tlx^p^ + Ca;a;a;'^""+V^""'^^ ■ (4) 

For nucleons, charge independence is respected by choosing a„„ = = a„p = a^N, with 
similar relations among the symmetry and repulsion coefficients. Thus, neutrons and protons 
can be treated as members of a single species with the appropriate symmetry term, and we 
can then retrieve the potential energy density of nucleons, similar to the form presented by 
LS (who denote the power in the repulsion term by 6): 

Epotin, z,p) = ^ a^^in + zfp"^ + b^Nin - zfp^ + CNN{n + zY^'^ p^'^^ . (5) 

LS use a slightly different expression, a'^^{n + z^p^ + h'^^nzp^ + CNN^n + zY~^^p^~^^, which 
requires a rearrangement of the coefficients in Eq. (^). In reaching Eq. (^, {n — z) repulsive 
contributions were dropped so that the short-range repulsion term depends only on the total 
nucleon density. 

We now turn to high density matter, where the hyperon species A, S and S are considered 
in addition to the nucleons. The full potential energy density expression will include several 
new terms, for the different baryon-baryon combinations: 

1 



potip) = ^ [aNNN^P^ + hNN{n - zfp^ + cjviviV'+V'^'J (6) 

A^ A 



+ aAHHAp^ + CAH (^A^S + a) p<+^ 

+ a^NNT.p' + 6Eiv(n-^)(E--E+)p2 + c^^(^j^L_n''T. + ^^A_s^iv)p^+i 
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+ 2 



where we define for each baryon species x — Px/ p, and introduce a shortened notation of 
N — n + z; S = S° + S~; E = E+ + E° + E~. There are no standard (isovector) symmetry terms 
for the A and E° baryons, since they have zero isospin; a nonstandard (isotensor) symmetry 
term for E's is disregarded. We denote the exponent of the density-dependent repulsion term 
for the nucleon-nucleon interaction by 5, following LS, while all other repulsive interactions 
are assumed to have a common exponent 7, although further diversity is, of course, possible. 

The kinetic energy density of the baryons contains the non-relativistic fermion kinetic ener- 
gies: 



£kin{{Xi},p,T) = J2 



n Ti(xi,p,T) 

2m* 



(7) 



where m* is the density dependent effective mass of the baryons of species i (with c — 1). The 
entire temperature dependence of the baryon thermodynamic functions is implicitly included 
in this term through the Fermi functions 



Ti{Xi,p, T) 



27r2 



2m*T 



5/2 



F3/2 {r]{xi,p,T)) 



(8) 



where the temperature is in units of energy (setting /c^ = 1). F3/2 is the corresponding Fermi 
integral, generally defined as 



(9) 



The parameter r]i is directly related to the temperature and density through the inverse 
relation 



77,(x,, p, T) = F-/^ 2^Xip {h^l2m*T) 



3/2' 



(10) 



In the case of zero temperature matter the Fermi function reduces to the simple expression 
Tj = I (Sn'^Xip)^^^ Xip. Typical nuclear models predict that in the depth of bulk matter the 
effective mass of baryons, m*, will be lower than the bare mass m^. The effective masses 
may be adequately represented in the Skyrme fashion 



1/m* = l/rrii (5ip 



(11) 



The scaling parameters /3j arise, however, from the model of the baryonic forces, and cannot 
be assumed independently. The effective equation has no means to fully relate consistent 
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effective masses and forces (unlike RMF models where effective masses result naturally), and 
this point must be considered separately. 



Finally, the contribution of the baryonic mass density is obviously given by 



(12) 



The chemical potentials of the baryons are also immediately related to the parameter rji : 



d{etot/p) 



dxi 



ViT + 



d{{epot + )/p) 



{n},p 



dxi 



(13) 



{ti},P 



In the special case of zero temperature matter rji is undefined, and the chemical potentials 
are found by analytically differentiating d{etot/ p)dxi. 

Lepton contributions to the energy density are calculated through the complete (relativis- 
tic) Fermi integrals. For sufficiently high temperatures (T > 1 MeV), the thermodynamic 
functions for electrons and neutrinos (when assumed to be trapped in the matter) can be ap- 
proximated by the ultra-relativistic limit. Keeping the lowest order mass corrections results 
in the following equations for the chemical potentials |1TT| : 



pYi 



9i ( Pi\ 



3 r 



TT 



-m, 



(14) 



where mi is the lepton rest mass (again in units of energy), gi is the spin degeneracy and 
Yi is the net lepton fraction per baryon (the number of particles minus the number of anti- 
particles). The chemical potential is then solved by 



pi = r -q/r ; 



(15) 



q^ + t^ + t 



nl/3 



t = M^h^pYi/gi ■ and q 



-vr 1 



-m, 



The contribution of the ultra-relativistic leptons to the total energy density for each lepton 
species I is given by: 



87r2 



vr^T^ / 7 

— TT^T^ 



Pf V15 



1 

-m 



(16) 



Obviusly, in zero temperature matter all formuals for the lepton thermodynamic functions 
reduce to simple zero-temperature fermion gas functions. 

We note that the muon cannot be treated as ultra-relativistic, and its thermodynamical 
functions must be calculated through the corresponding Fermi integrals. However, in the 
case of neutrino-trapped matter mouns are almost extinct (unlike the neutrino-free matter 
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case, where muons form at significant fractions), and we shall disregard them for the case of 
finite temperature matter. 

The construction of the equation of state is completed by the calculation of the various 
contributions to the pressure in the matter. The pressure is derived from the total energy 
density by: 

P = Y1 ^iP^^i + H ^lP^^l ~ f = Y1 ^iPPi + ^iPPi + Ts-e , (17) 

i I i I 

where / is the free energy density, f = e — Ts, and s is the entropy density. Again, the 
pressure is composed of a baryon contribution and a lepton contribution, easy to derive from 
Eqs. (|6D- (p!3D and (p!^)- (|T6|) respectively. In finite temprature matter the entropy per baryon 
(the commonly used variable in hydrodynamical simulations), s, for each baryon species is 



and for each of ultra-relativistic lepton species is 

' „ „ 1 _\i 

(19) 



. _ 1 giTfii 



2 / 7 1 

V15 2 ' 



Finite temperature calculations must include a photon gas /i^ = contribution, with: 

_ ^^^^ p _ 1 A ~ - ^^^^ 

^7 ~ ^g^3 ; Pi — 2^7 ^"^^ ^1 ~ ~Ij~J^ ' 



3 Equilibrium Conditions 



The composition of the high density matter is obviously constrained by charge neutrality 
and baryon number conservation. Further constraints are imposed by the assumption of 
complete beta-equilibrium as the time scales of the weak interactions are generally much 
shorter than dynamical time scales in the relevant stellar systems. It is this constraint that 
sets most of the equilibrium conditions, as it forces all processes to take place in absolute 
chemical equilibrium. 

The balance between baryons and leptons is maintained through beta decays. In nuclear 
matter, this is simply the equilibrium of the process p + e~ ^ n + Ug. This equilibrium 
establishes the relation 

fJ'p + fJ'e = fl'n + f^u, ■ (21) 

Strange baryons may form in the matter through various types of weak interactions. For 
example, the A hyperon can appear both in leptonic processes, p + e~ — > A + z/g, and in 
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baryonic processes such asn + n^A + ra (the inverse of which occurs in hypernuclei) . There 
are also other possible processes, involving pions and kaons; however in this study we assume 
that meson condensates do not form in the matter (see discussion later in this section) and 
are therefore excluded. Note that the balance among strange baryons is maintained mainly 
through strong interactions like A + n — + P, or A + A + p. The resulting 

equilibrium conditions are simply that the chemical potentials of all baryons depend only on 
their electric charges (that determine the possible combinations in the various reactions): 

/Un = /LtA = fij^o = fi^o ; fip = ; Ats- = f^s- = '^f^n - f^p ■ (22) 

It is noteworthy that the chemical potentials of the negatively charged hyperons are condi- 
tioned by a larger value than the neutral ones (since dense matter is typically dominated 
by neutrons). For example, studies of cold high density matter usually find that the first 
hyperon to appear is the S~, closely preceding the lower- mass A due to the charge difference. 

The last fermion players in the high density matter are the muon and the muon neutrino. 
When the electron chemical potential exceeds the muon mass, muons and muon neutrinos 
can be formed by e~ — > /x" + z/g + ^^t, so that the equilibrium condition for the muon chemical 
potential is: 

Eqs. (pTD-(p3D impose eight conditions on twelve variables, so that another four conditions 
are needed to uniquely determine the complete equilibrium composition of matter at given 
density and temperature. Two further conditions arise from charge neutrality and baryon 
number conservation. These conditions constrain the species fractions {xi} and {Yj} by 

Y,Xiqi + J2Y,q,=0 (24) 

i j 

where {qi} are the electric charges, and 

E^* = l- (25) 

i 

Finally, two more conditions result from the physical assumptions about the neutrino diffu- 
sion time scale. Core-collapse and proto-neutron star matter may be assumed to be opaque 
to neutrinos, so that neutrinos are trapped and the total lepton fraction must be kept con- 
stant. As the trapping in supernova occurs when the collapsing core reaches densities of 
about 10^^ gm/cm^, where no muons exist, the two resulting conditions are: 

Ye + Yu^ = Yi = const ; Y^ + Y,^ = (26) 

where Yi is determined at the onset of trapping at about 0.35-0.4, with Y^/Y^^ ^5 — 7, 
depending on the exact trapping density (see Bethe's detailed review [jl8|). As metioned 
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earlier, we will ignore the minute muon production in proto-neutron-star in this study. Neu- 
tron star matter, on the other hand, is completely transparent to neutrinos (the transparency 
settles after about 50 seconds of proto-neutron-star evolution [0), and the neutrino chem- 
ical potentials are set to zero. Muons are then easier to create (the threshold condition for 
muon creation is simply fie = iTif^), but the total lepton density is not constrained. Neutron 
star matter is considerably deleptonized with respect to the proto-neutron star. The two 
equilibrium conditions are: 

fJ^e = (J^fi ; fJ'u, = fJ'u^=0 , (27) 

and all the chemical potentials can be expressed as linear combinations of the electron and 
neutron chemical potentials (associated with the electric and baryonic charges, respectively). 

Solving the complete set of conditions yields the equilibrium fractions of the baryon and 
lepton species for any given density and temperature. The effective form of our EoS enables 
a low-order Newton- Raphson iteration process, which is quick and efficient. The equilibrium 
compositions and the resulting thermodynamic properties are then calculated for the large 
variety of strong interaction parameterizations discussed below. 

Other particle candidates for appearance in the high density matter are the A isobars and 
TT and K mesons as condensates. Current theories predict ttN S-wave repuslive forces which 
exclude pion condensates, and indirectly, make the formation of A's unfavorable as well p, 
T0[] . However, as shown by Weber and Weigel [^], in Relativistic-Hartree-Fock calculations 



some evaluations of the pNA coupling enhance the formation of the A~. Kaon and anti- 
kaon condensation has attracted considerable attention [^, but some recent investigations 
(see, for example, Refs. ^) do not find kaon effective masses dropping low enough at 
relevant densities. Clearly, the formation of hyperons in the matter will further suppress 
meson condensation p, ^. We therefore opt not to include A's and meson condensates in 
this study; further generaliztion is, of course, possible, where the model dependence of the 
above arguments may be examined. A final issue is that of the transition into a deconfined 
quark phase: at some high density the hadronic identities of the quark states must begin to 
fade, as the baryons will overlap. In this study we do not treat quark degrees of freedom, 
and allow the baryonic effective equation of state to dominate up to Sps- Quark degrees of 
freedom, however, might come into effect at lower values of density (an especially exotic case 
is the possibility of stable strange quark matter, which leads to the existence of strange stars 
[BH]). The issue of quark deconfinement is defered to a future work. 
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4 Choice of Parameters 



The local potential approach requires four parameters for every combination of baryon- 
baryon interactions in Eq. (^: the coefficients a,b,c and the density power dependence 7 
{6 for the nucleon-nucleon case). The effective masses of the baryons, or rather their density 
dependence through the parameters Pi defined in Eq. (|TT|), are required as well. 

In principle, all parameters of the local potentials could be determined experimentally, but in 
practice there are large uncertainties due to limited experimental data, and some values must 
be drawn from theory. The experimental and theoretical data used as basis for evaluating 
the parameters are discussed briefiy in the proceeding paragraphs. 

4.1 Nucleon-Nucleon Interactions 

The properties of the nucleon interactions at the nuclear saturation density, ps, are measured 
to a reasonable precision in finite nuclei. The experimental value of this density is estimated 
to be in the range 0.145-0.17 fm~'^, and in this study we follow LS by employing the Skyrme 
r force value of ps = 0.155 fm~^. 

As shown in LS, the values of the parameters otvat, fcATAr, ctvat and S are evaluated from 
the binding energy, symmetry energy and incompressibility of saturated symmetric nuclear 
matter, derived from typical experimental results and the Skyrme model (see Section 2.2 
of LS for detail). However, the parameters aj^NjCNN' and 6 strongly depend on the nuclear 
incompressibility which is poorly determined in finite nuclei, and is difficult to match with 
theoretical predictions (which generally yield too high incompressibility). Thus, quantifying 
the properties of the nucleon-nucleon interactions at high densities leaves a large range of 
uncertainty, especially for 6, which as a result may be set practically arbitrarily between 
unity and about 2. For example, LS use the Skyrme value of the incompressibility Kg = 375 
MeV, leading to 6 = 2.002, but in a later work values down to Kg = 180 MeV, resulting 
in 6 = 1.07 are also considered. In order to investigate the model dependence due to this 
parameter, we will examine the values 5 = 2, | and |. The appropriate values for the 
parameters a^^jb^N and for these values of 6 are presented in Table 1. 

4.2 Hyperon-Nucleon Interactions 

Experimental data of effective nucleon-hyperon interactions are confined to hypernuclei (for 
a recent review, see [^). This clearly limits the interpretation of the data, and requires a 
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combination of experimental and theoretical analysis. 

There is a considerable body of information on binding energies, -Ba, of A's bound in various 
single particle orbitals in hypernuclei with total baryonic number from A = 3 to A = 208. 
By analyzing these data in a Skyrme-Hartree-Fock picture MDG determined a well depth of 

- Va{pn = Ps) - 27- 28 MeV , (28) 

defined as the potential seen by a single A in nuclear matter at saturation density, which 
is roughly | of the nucleon well depth in a similar parameterization. Constraining the 
local potential by the well depth at nuclear saturation density is the equivalent of using the 
bulk binding energy to constrain the nucleon-nucleon interaction (related to the well depth 
through the Fermi kinetic energy). By combining the data for many hypernuclei, MDG used 
the A dependence to evaluate the A-nuclear potential for densities p < Ps and fitted a local 
potential of the form Va{pn) = clPn + cpj^ to the experimental results, deducing a and c for 
a given 7. MDG achieved reasonable fits for the single A orbitals in ordinary hypernuclei, 
but - again - could not constrain the shape of the potential for the high density regime, and 
therefore could not uniquely determine the exponent 7 of the repulsion term. Several possible 
choices for the parameters a, c and 7 are presented in Table 2 below. Figure 1, reproduced 
from MDG, demonstrates the character of the uncertainties of the baryon interactions, as 
inherent in the local potential formulation. The single A local potential in nuclear matter, 
V\{pn), is plotted against the nuclear matter density, for different choices of the exponent 7. 
All the potentials are constrained by a well depth of about —28 MeV at nuclear saturation 
density and the A-hypernuclear binding energies for lower densities, so that all the plots 
yield similar results for p < Ps- However, the unconstrained behavior at high density is quite 
divergent, depending strongly on the value of 7. For example, the cross-over density, pco 
for which Va(Pco) = varies from about pco = '^Ps for 7 = 2 to pco = 3.5ps for 7 = |- It is 
noteworthy that this local potential does reproduce the calculated single particle potentials 
found in RMF models l3( 



Experimental data of S hypernuclei are more ambiguous, partly because of the strong SN— 
AN decay channel. It is common to assume that the E well depth is equal to that of the 
A ||2^. The isospin dependent component of the interaction of a S hyperon in nuclear 



matter (S-nm) is usualy evaluated from theory, typicaly predicting a E isovector coupling 
constant (due to the p meson field in RMF models) of the same order of magnitude as for 
nucleons. Possible choices for the parameters a, 6, c and 7 then follow from Tables 1 and 2, 
and are presented in Table 3. However, a recent analysis |]13|, 0] of theoretical work and 



experimental data concerning S~-atoms suggests that there is a strong isoscalar repulsion 
in the S-nm interaction. These studies, which also determine the isovector component of 
the E~-nucleus interaction, indicate that E baryons are strongly repulsed in the nuclear 
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interior. A corresponding density dependent local potential has also been proposed (again 
with difficulty in determining the high power density dependence) for S~ hyperons: 



2/iV^E(p 



N 



An 1 + 



m. 



ho + B, 




+ 




{n-z)\pN , (29) 



where /i is the S-nucleus reduced mass (for nuclear matter it is simply the S mass) and m 
is the nucleon mass. The (60, -Bo) term describes the isoscalar component of the interaction, 
and with Bq <Q and |-Bo| > ^0 a strong repulsion is produced in the nuclear interior. The 
second term introduces isospin dependence, reversing sign for the S"*". The real part of this 
potential may be taken as the S-nm interaction for equilibrium purposes, as the imaginary 
part represents mainly the SN^ AN channel. A representative set of parameters for Eq. 
( PUj) is given in Table 4. We investigate both of the possibilities listed in Tables 3 and 4 for 
the S-nm interaction in constructing the equation of state. 

In emulsion experiments with K~ beams there are a few events attributed to the formation 
of S-hypernuclei. Although the evidence is far from being compelling, the data can be 
interpreted consistently in terms of a potential well depth of 

- ■^h(P7v = P.) ~ 20 - 25 MeV . (30) 



Such a result has also been derived in one-boson-exchange (OBE) models [^]. The isospin 
component of the SN interaction is very model dependent, ranging from negligible |2^ (an 
order of magnitude smaller than that observed for nucleons) to equal to that of nucleons 
[||. Possible choices for the parameters a, c and 7, assuming h~_N = are listed in Table 5. 
In deriving these parameters we assumed —V~,{pn = Ps) = 24 MeV for symmetric nuclear 
matter (where symmetry terms vanish) and that for a given value of 7 the crossover density, 
Pco, for which Vs{pn = Pco) = is equal to that for A in nuclear matter, thus simulating 
universality in the high density regime. 



4.3 Hyperon-Hyperon Interactions 

The experimental data regarding hyperon-hyperon interactions is extremely scarce. There 
are a few events which may be interpreted by the creation of AA hypernuclei (see Ref. 
p5|]). Analysis of these events requires a rather strong AA attraction, which is not easily 



reproduced by theory . Recent studies with RMF models remedied the hyperon-hyperon 
interaction by adding two mesons ^ which couple only to the hyperons. Schaffner et al. 



23] constructed single particle potentials on the basis of OBE calculations of the Nijmegen 



group, and estimated the well depths of hyperons in hyperon matter as follows: 

-Vh,iph,=Ps) = ^OMeV , (31) 



13 



implying a universal hyperon-hyperon interaction. Similar to the case of the S-nm interac- 
tion, the coefficients flfe^/ij (^hih2 can be determined for a given value of 7 by requiring that 
for symmetric matter the potential depth at a density equal to the nuclear saturation den- 
sity be equal to the value assumed in Eq. (|3T|), and that the crossover density for which the 
potential is zero, be equal to that of the A in nuclear matter (A-nm) case. The appropriate 
values for the inter- hyperon interactions are given in Table 6. 

Generally, the symmetry coefficients hxy in Eq. (H) are independent of 7, and must be de- 
termined separately. The isospin dependent components of the hyperon-hyperon interaction 
must be drawn from theory. As mentioned earlier, the isovector coupling constant of the 
S baryons is typically assumed to be up to twice that of the nucleons, while that of the S 
baryons is unclear. We assumed that only the N-nm, S-nm and S-Em interactions have a 
nonzero isospin component, and the corrseponding values of the symmetry coefficients are 
presented in Tables 1, 3 and 6. 



We conclude this discussion by reemphasizing the large uncertainties surrounding the baryon- 
baryon interactions, especially regarding the density dependence. Constraining the local 
potentials in nuclear matter by assuming a given well depth at saturation density does not 
limit the shape of the potential at higher densities, which is handled in the local potential 
forms of Eq. (P) by the parameters 7 and 5. We therefore treat the exponents 7 and 5 as 
free variables, and simply assume a universal nature to all multibody forces, so that we may 
also set 7 = 5. 



The evaluation of the baryon effective masses is a fundamental problem in any analysis of 
supernuclear density matter. The values of effective masses are needed for the kinetic and 
mass energy density terms of Eqs. ^ and (12) respectively, and nuclear matter calculations 
typically predict [|16| that the nucleon mass in bulk matter is smaller than the bare mass 
by a factor of m* /m ^ 0.7 at saturation density, and will continue to decrease for higher 
densities, due to both interaction and relativistic effects. The analysis of the effective masses 
of all the baryons in yields similar results for nucleons and hyperons alike. 

The local potential has no mechanism to generate both the interactions and the effective 
masses consistently. LS chose to set m*/m = 1, which is a reasonable approximation for 
densities up to about twice the nuclear saturation density (and may also be justified by some 
uncertainties regarding effective masses at the Fermi surface [^). At higher densities such 
an approximation becomes questionable. While the explicitly density dependent interactions 
in the potential model may compensate for some of the expected decrease in the effective 
masses (LS include a further explicit dependence of the forces on the assumed values for the 
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effective masses - see Eqs. (2.20-21) in LS), relativistic effects are not taken into account. 

For lack of a better model, in this work we will use m* /m = 1. An obvious flaw in setting all 
effective masses equal to the bare ones is that mass differences between the various types of 
baryons remain constant. Theory clearly suggests otherwise (see, for example. Figs. 4 and 
7 in P]), and changing the mass differences will clearly affect the equilibrium compositions. 
For example, theory predicts that the effective mass of nucleons in nuclear matter at p = Ps is 
^Ar/'"^Ar~ 0.55 — 0.7, while that of the A decreases only to m\/m\^0.8 fl^; hence the mass 
difference will be larger than for the bare masses (and A's will form in nuclear matter at 
higher densities than for m*/m=l). The sensitivity of the results to this crude assumption 
can be examined by using nonzero positive values for the coeffecients (3i of Eq. ([TT|) , which 
will yield effective mass functions similar to those found in RMF models (for example, in 
Refs. 1,1). 



5 Results: Equilibrium Compositions and Equations 
of State 

The basic features of the dense matter are the different particle fractions and the equation of 
state, both of which are dependent on the baryon density of the matter, and of the modeling 
of the strong intearctions. In this Section the equilibrium compositions and equations of 
state are examined and compared for the different cases of candidate hyperon species (no 
hyperons, A's, S's and H's, and only A's and H's but no S's when assuming a repulsive E-nm 
interaction). The sensitivity to the choice of parameters of the baryon-baryon interactions 
is also examined, by setting different values for the exponents 6 and 7. 

We shall focus on the simpler case of neutron star matter, which can be well approximated 
as being at zero temperature, where all thermodynamic functions are easily found. Some 
results regarding proto-neutron-star (core- collapse) matter are presented in the latter part 
of this Section, where the assumed thermodynamic conditions are detailed. 



5.1 Neutron Star Matter 

Neutron star matter may be well approximated as having zero temperature, since evolved 
neutron stars have internal temperatures considerably below 1 MeV. The internal energy 
density e and the free energy density / are identically the same, and the kinetic energy 
terms in Eqs. (0,|) are all simply related to the Fermi momenta. 
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Shown in Fig. 2 are the particle fractions versus the baryon density for matter containing 
neutrons, protons electrons and muons, calculated by using the parameters for 5 = 2 (see 
Table 1), and ignoring the possible appearance of hyperons. Compositions are calculated 
up to about 8ps, where quark degrees of freedom are assumed to come into effect, and the 
validity of a baryonic model of the matter becomes questionable. These plots are very similar 
to those calculated in various models of nuclear matter for neutron stars. The main features 
of such matter are: 

1. Neutrons are by far the dominant species in the matter. 

2. Both electrons and muons are present in the matter. At increasing densities the chem- 
ical potential of the leptons rises (due to the rising Fermi energies), and the muon 
fraction approaches that of electrons, as the Fermi energies drive the chemical poten- 
tial up so that the mass difference is of lesser significance. 

3. The balance of nucleons and leptons is determined by the difference /Xn ~/^p (commonly 
denoted as fi) and by /ig. Equilibrium demands that both be equal, but since /i„ — /ip 
has a steeper density dependence, there is some decrease in the neutron fraction and 
some increase in the proton and total lepton fractions as a function of the total baryon 
densities. 

These features are common to a variety of nuclear matter models and are only weakly model 
dependent. This results in similar graphs for the other parameterizations; although specific 
details are, of course, model dependent (for example, the exact proton fraction which is 
highly important for direct Urea cooling rates). 

The picture clearly changes when hyperons are taken into account. Figures 3a, 3b and 3c 
show the various particle fractions as a function of the total baryon density for 7 = 5 = 2, | 
and I respectively. In all cases the hyperon-hyperon well depth was set as —40 MeV, and the 
E-nm interaction is assumed to be equal to the A-nm one (except for the symmetry term, 
which is nonzero for E's and nucleons). The appearance of hyperons affects all the main 
features of the equilibrium compositions. The key effects to be noted are: 

1. In all cases hyperons appear at about twice the nuclear saturation density, and at den- 
sities of about 2.5ps they sustain a significant fraction of the total baryon population. 

2. The first hyperon species to appear is the S~, closely followed by the A. The negative 
charge of the S~ outweighs the 80 MeV mass difference, as a result of the more lenient 
condition, Eq. (P^D, for its chemical potential. However, the formation of hyperons 
is quickly moderated by the charge dependent forces, which disfavor an excess of S~'s 
over S"'"'s, and a joint excess of S^'s and neutrons. Thus the S~ fraction saturates at 
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about 0.1, while the A's, free of charge-dependent forces, continue to accumulate until 
the repulsive (multibody) forces bring saturation on them as well. 

3. All other hyperon species follow, and appear one by one in the matter. With the 
assumptions mentioned above, E's generally appear before the S's due to the large 
mass difference, but the S~ becomes favored due to its negative electric charge, and 
quickly becomes abundant in the matter. Since we assume no charge-dependent forces 
for the S baryons, the E" fraction does not saturate quickly, and its fraction overtakes 
that of the E~ at high densities. 

4. A unique aspect of hyperon appearance in all cases is the immediate deleptonization 
of the matter. Leptons are rather expensive in terms of energy density (and pres- 
sure), surviving in nuclear matter only because of the need for charge neutrality and 
because the charge-dependent forces impose a significant proton fraction. Hyperons 
create an option for lowering the neutron excess free of lepton formation, and the neg- 
atively charged hyperons allow charge neutrality to be maintained within the baryon 
community. Thus, the increase of the lepton abundance is cut short by the formation 
of the hyperons and, in some cases, the appearance of the S~ is followed by a very 
powerful deleptonization. For 7 = (5 > | the muon population is even completely ex- 
tinguished (obviously, the E~s cannot achieve this alone because of their saturation 
fraction, imposed by the symmetry forces) . The case of 7 = 5 = | is somewhat eccen- 
tric, maintaining a high lepton content along with the rather moderated appearance of 
hyperons. This is a consequence of setting ^ = |, which represents a too mild increase 
in the nucleon-nucleon repulsion with increasing density, thus suppressing hyperon pro- 
duction. Such a composition seems to be unlikely, since a combination of hyperons and 
a high lepton content will create very rapid cooling rates [32-35] . 

Some model dependence docs exist, of course, as weaker nucleon repulsion is necessarily 
accompanied by less hyperon production. For S — 2, hyperon production is so favorable that 
at high densities A's, not neutrons, are the most abundant species, while for 5 = | nucleons 
dominate the matter even at 8ps- The appearance of hyperons other than E~ and A is also 
somewhat model dependent: the more powerful the high density repulsion (among nucleons 
and hyperons) , the formation of E*^ and even E+ becomes more favorable at sufficiently low 
densities, where the more massive S~ is still absent; for weaker repulsion, it is the E~ which 
appears first of the three, at densities where the effect of its negative charge is sufficient 
to overcome the 120 MeV mass difference. S°'s appear only when repulsive forces become 
dominant again and the introduction of more species is thus more effective in lowering the 
potential energy density. 

The reason for the hyperon formation density being only weakly model dependent is that 
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when hyperons are excluded, the nucleon chemical potentials in nuclear matter rise steeply 
with density. All models predict that at about 2ps the nucleon-nucleon interaction is far 
beyond maximum attraction, and in some models it is already repulsive, so that lowering 
the partial density of the nucleons is clearly favorable. Thus, although the exact density for 
hyperon appearance is model dependent, it does not vary by more than a few hundredths of 
a fm~^. We note that this is also the reason for the quick rise in the hyperon fractions with 
density after formation. The only balance to this trend are the kinetic energies, which are 
far outweighed until the hyperon fraction rises to the order of about 0.1. This behavior is 
less pronounced for the E's than for the A's (and in our model, also for the S's) because of 
the nonzero symmetry term, which disfavors a fast build-up of one species. 

When using the effective equation of state, the effects discussed above tend to overshoot, 
creating non-monotonic irregularities in the curves (most notable in the proton and lepton 
fractions). This (unphysical) effect is moderated to large extent in RMF models, which add 
further control through the effective masses which are density and species-fraction dependent. 

Generally, Figs. 3a and 3b show a fair agreement with the results of the equilibrium com- 
positions found in the mean field calculations of [3-5]. This agreement suggests that the 
effective equation of state can indeed be considered as a trustworthy tool in the analysis 
of high density baryonic matter. The somewhat eccentric case of 7 = 5 = | examined in 
Fig. 3c, however, is qualitatively different from the RMF results. It seems that this case 
should be treated with some suspicion, although we will keep it in the analysis below. It 
is noteworthy that the quantitative details of the compositions are model dependent in the 
RMF calculations as well. For example, the deleptonization rate and other features of the 
PL-Z model in Ref. resemble our 7 = 5 = 2 case, while those of case 2 in Ref. ^ (which 
is relevant for our discussion) bear higher resemblance to the 1 = ^= \ case. 

We now turn to the case of strong isoscalar repulsion in the S-nm interaction, suggested in 



1^, Using the local potential of Eq. ( pOj) with the parameters listed in Table 4 (or any 



of the other sets of parameters in that fit the experimental data of S~-atoms) yields 
a repulsive potential of the order of tens of MeV's for E's in nuclear matter at saturation 
density, and even greater repulsion at higher densities. S hyperons are thus disfavored in 
nuclear matter, and the obvious byproduct is the exclusion of these species in the high 
density matter under discussion. 

Figures 4a, 4b and 4c show the various particle fractions as a function of the total baryon 
density for 7 = 5 = 2, | and | respectively, when the hyperon-hyperon potential well depth is 
set as —40 MeV, and the S-nm interaction is constarined by fitting to the E^-atoms data. 
Clearly, S hyperons are absent, but other than this, the key features of hyperon formation 
are unchanged: 
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1. Hyperons appear again at about twice the saturation density, and at about 2.5ps the 
hyperon fraction is already significant. In this case, however, the entire effect is due 
to formation of A's. This quantitative result is practically model independent. The 
appearance of the A's halts the increase in the lepton population, and when S~'s are 
added further deleptoniztion occurs. The A actually appears at lower densities than it 
did when E's were present, since now it has no competition from the E~. 

2. The S~ appears as well in all cases, at shghtly lower densities than for the corresponding 
case when competing with E's. The formation of S° is model dependent, again, since 
stronger high density repulsion favors the introduction of more species. 

The total hyperon fraction is again highly model dependent as are the specific details of the 
deleptonization. 

The sensitivity of the results to the assumed value of the hyperon-hyperon well depth, is 
examined by setting its value to zero - "turning off" the hyperon-hyperon interaction (with 
the exception of the E-Em symmetry term, which is kept unchanged). Fig. 5 presents the 
equilibrium compositions found for the ^ — 5 — 2 case when all baryon species are allowed 
(assuming the E-nm interaction to be equal to the A-nm one) for a zero hyperon-hyperon 
well depth (while all other parameters are unchanged). The considerable similarity to the 
compositions for the corresponding case with a relatively large hyperon-hyperon well depth 
(Fig. 3a) is clear. There is a shght moderation in the rate of accumulation of hyperons 
(missing the attractive interactions among hyperons of the same species at small partial 
densities) while the appearance of the more massive hyperons is somewhat delayed, because 
there is also no repulsion among the abundant hyperons. These effects are visible but 
clearly are of second order, since the matter is dominated by nucleons in most of the density 
range, so it is the nucleon-nucleon and hyperon-nucleon interactions which determine the 
equilibrium compositions. The sensitivity to the hyperon-hyperon interaction is even weaker 
for lower values of 5 and 7, as the nucleon fraction in these cases is larger. A similar relative 
insensitivity to the hyperon-hyperon well depth is found also when E's are excluded from 
the matter. 

The solution for the equilibrium compositions enables the determination of the thermody- 
namic functions, explicitly the energy density and the pressure, for constructing the EoS. 
The key point regarding the effect of introducing extra particle species into the matter in this 
respect is the softening of the EoS. This softening is the consequence of turning potential 
and kinetic energy of the already present baryons into masses of the new species, and in 
reducing repulsion (and sometimes even gaining attraction) in the potential energy terms, 
due to the lowering of the partial densities of the abundant species. These two effects lower 
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the energy density (so that the mass-energy density at a given baryon density is somewhat 
less than the corresponding energy density of nuclear matter at the same baryon density), 
but the main effect involves the pressure, since storing energy in mass is cheapest in terms 
of pressure. The decrease of the lepton population does add further softening, but this effect 
is of second order. 

The nature of the EoS is thus dependent on the number of species, in addition to the 
obvious dependence on the baryon-baryon interaction model. It is instructive to begin with 
the resulting EoS for nuclear matter for the three cases corresponding to 5 = 2, | and |, 
plotted in Fig. 6. Since no hyperons are considered, the resulting EoS is a direct function of 
the nuclear interaction model, namely, of the value of 6. The differences between the three 
equations are clearly visible. It is noteworthy that the 5 = 2, | equations break the causality 
limit {{dP/de) < c^) at some density, where the plots are made dashed. This is a known 
flaw in effective models for the nuclear interactions, sometimes amended by introducing a 
finite-range into the local potentials. 

The properties of the EoS change qualitatively when hyperons are taken into account: Fig. 
7 shows the resulting EoS for three types of matter with 5 = 7 = |, for a hyperon well depth 
of —Vh^h^ = 40 MeV. The three cases are 1) nuclear matter (n-p-e-yu); 2) baryon matter 
including A's and S's but no S's (assuming a repulsive S-nm interaction [0); and 3) all 



eight species of baryons appear. The three plots demonstrate clearly the softening of the 
EoS when hyperons are included, this softening becoming more pronounced when the S's 
are also included. 

We must note that the numerical peaks in the particle fractions discussed above are unfortu- 
nately carried over to the EoS, when the unphysical fast rise in the A fraction (and sometimes 
the S~ fraction as well) yields a slight dip in the pressure, so that over a small region the un- 
physical result of dP/ de < is obtained. By employing a mathematical smoothing scheme, 
interpolating the pressure- density dependence in these small regions, we are able to amend 
this shortcoming, particuUary for the case of 7 = 2. 

Fig. 8 compares the resulting equation of state for three cases of the baryon-baryon inter- 
actions, i.e. 7 = 5 = 2, 1 and |, and in all cases the hyperon-hyperon well depth is assumed 
to be —Vhih2 =40 MeV, and the E-nm interaction is assumed to be equal to the A-nm one. 
It is noteworthy that all three equations of state are quantitatively similar, indicating a 
rather weak dependence on the precise description of the short-range forces represented by 
the exponents 5 and 7. A similar trend is found when a highly repulsive S-nm interaction 
is assumed, with all equations being somewhat stiffer than the coresponding cases where 
the S's are allowed in the matter. Thus, an important feature of our results is that the 
equation of state is much less sensitive to the details of the baryon-baryon interaction, with 
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the obvious exception of the two extreme possibihties for the E-nm interaction. 

The nucleon-nucleon interaction, which determines the stiffness of the nuclear-matter equa- 
tion of state, also determines the rate at which the nucleon chemical potential rises with 
increasing density. A greater rise in the nucleon chemical potential is the incentive to en- 
hanced formation of hyperons in the matter, which in turn softens the total equation of 
state. Our key observation is that the macroscopic properties of the equation of state are 
thus strongly coupled to the microscopic equilibrium compositions, and once hyperons are 
available the resulting equation of state is constrained to a significantly narrower range of 
values. We may conclude that hyperons serve as a "pressure control" mechanism in high 
density matter, and its nature depends mainly on the number of hyperon species available. 

Explicitly, as long as attractive nucleon-nucleon and nucleon-hyperon well depths are as- 
sumed, the resulting EoS is only mildly sensitive to the choice of 5 and 7. This is mainly 
reflected in the width of the "semi-plateau" region, where the pressure grows only very 
mildly: the greater the value of the exponents, the higher the relative gain \dV/dp\ in attrac- 
tion at low densities and the wider the density range where the hyperon appearance, mainly 
of the A, creates the " semi- plateau" region. On the other hand, if a large value is chosen for 
the exponents, there is also a steeper rise of repulsion at higher densities, which builds up 
the pressure faster. These two effects are balanced by the appearance of new species in the 
matter, when such species are available. New species will appear when the pressure (and 
internal energy) rise steeply with the baryon density, but their appearance is delayed when 
the species already present moderate the pressure build up, as happens in the " semi- plateau" 
regions. 

The softening of the EoS may be demonstrated by calculating the adiabatic index F = 
e/ Pd{P)/d{e). Fig. 9 shows the adiabatic index as a function of the mass energy density, 
for nuclear matter and for matter containing nucleons and hyperons, with 5 = 7 = |. The 
appearance of the hyperons is easily identified, as each of them lowered the adiabatic index 
as long as its fraction did not saturate. The A's and S's cause greater "dips" in the value of 
F than the E's, since their formation is not moderated by charge dependent forces. 

The sensitivity of the "pressure control" to the hyperon-hyperon interaction is found to be 
especially weak, since the matter is nucleon dominated in at least most of the relevant density 
range (and for 7 = 5 = | even for the entire range). The deeper Vhih2 is, there will be a stronger 
hyperon-hyperon attraction at low partial densities (further favoring hyperon production), 
but also a stronger hyperon-hyperon repulsion at high partial densities. Hence, through 
most of the density range the hyperon fractions quickly saturate so that their contribution 
to the energy density (and pressure) are relatively independent of the details of the hyperon- 
hyperon interaction. 
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We comment briefly on the issue of tlie effective masses. Several equilibrium compositions 
were recalculated with baryon effective masses m*/m < 1 derived by using nonzero values 
for the coeffecients /?j in Eq. (]TT|). We find that as long as the mass differences between 
the various baryon species do not change considerably from the case of bare masses, the 
equilibrium compositions are very similar to those found by using m*/m= 1. This is a natural 
result, since the chemical potentials are highly density dependent through the potential 
energy terms, and so the densities at which hyperons appear do not vary by more than 
a few hundredths of fm~^. The corresponding equations of state are more sensitive to 
the values of the effective masses, requiring some careful consideration. We find that the 
pressure vs. baryon density is only mildly dependent on the precise choices for m*/m, 
and the basic feature of hyperon induced "pressure control" is maintained. Since the mass 
density gives a large contribution to the energy density, the pressure vs. energy density 
relation is more complicated when effective masses are used (and, indeed, LS scale their 
energy density relative to the neutron effective mass). However, if the energy density is 
scaled so that for nuclear matter at nuclear saturation density it yields the typical value of 
£{Pn = Ps) ~ 2.4 X 10^'' gm/cm'^, the resulting equations of state are once again approximately 
similar to those found for m* /m = 1. In summary, we believe that the equations of state 
found using m*/m = 1 reasonably reflect the gross properties of the equations of state that 
are calculated with finer approximations regarding the effective masses. 



5.2 Proto-Neutron Star Matter 

The physical conditions in proto-neutron star matter basically differ from those in fully- 
evolved neutron star matter in two main features: the temperature is of the order of 10 MeV 
and more, so it must be treated as finite, and - to a reasonable approximation - the matter 
is opaque to neutrinos, trapping them and keeping the total lepton fraction, Yi, constant. 

Both of these general considerations affect the equilibrium composition as well as the equation 
of state of the high density matter. Finite temperature means also finite entropy, which will 
stiffen the EoS of matter for a given composition |^ and broaden the population levels of the 
baryon species, thus lowering the chemical potentials. A constant lepton fraction imposes 
a much higher electron fraction than in neutron star matter, clearly affecting the baryon 
species composition, and through it the EoS. Generally, studies of core-collapse and proto- 
neutron matter with finite temperature and a constant lepton fraction point to the latter as 
the more important of the two. Its value is determined by the onset of neutrino trapping in 



the collapsing core of the supernova, and is usually found to be YJ = 0.3 — 0.4 [18|; in this 



Section we use the value of 0.36. The temperature is typically constrained by the entropy 
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per baryon, s, which tends to be constant throughout the star. We use a temperature of 
T = 20 MeV, which is about equivalent to s ~ 2 p. 

Once more we begin with nuclear matter, employing the different assumptions on the nuclear 
EoS. Fig. 10a shows the equilibrium composition of n-p-e-i/g, with the exponent S of the 
nuclear EoS set at (5= |- The basic feature of the equilibrium composition is a large electron 
fraction, imposed by the condition that Yg + Yu^ = ^i- The high electron fraction imposes 
an equal proton fraction, so proto-neutron star matter is much more symmetric at p ^ 
Ps {Z/ A 0.3) than neutron star matter {Z/A 0.1). It is also noteworthy that the neutrino 
fraction decreases with density (and, correspondingly, the electron fraction increases). This 
is the result of the equilibrium condition (|2l|), and since /i„ — /ip rises more rapidly than 
/^e — Pvei the equilibrium is maintained by lowering the neutron fraction. These trends are 
highly model independent, and similar equilibrium compositions are found for the other 
values of the exponent 5. 

Equilibrium compositions of finite temperature matter will always contain some finite pop- 
ulation of all particle species, so minute hyperon fractions exist even at ps- However, these 
hyperon fractions will build up to significant values only when the potential and kinetic 
energies of the nucleons make the transition of nucleons to hyperons energetically favorable. 
Hence, the general features of proto-neutron star matter will be qualitatively similar to those 
of neutron-star matter discussed above, and significant hyperon formation is expected. 

An example for the effects of neutrino trapping and finite temperature on the equilibrium 
compositions is plotted in Fig. 10b, for the case 7 = 5= | and —Vh^h2 =40 MeV, for Yi = 0.36 
and a temperature of T = 20 MeV. These compositions appear to be generally similar to those 
of the corresponding case for neutron star matter (Fig. 3b). Some noticeable differences do 
still exist, of course: the greater symmetry of the nuclear matter at low densities, imposed by 
the large electron fraction, reduces the neutron chemical potential with respect to neutron 
star matter at the same densities. As a result, significant hyperon formation is somewhat 
delayed in the neutrino trapped matter. The finite temperature slightly enhances this effect, 
since the nucleon chemical potentials drop below the corresponding value in cold matter. 
Details relevant to the hyperon formation are as follows: 

1. The formation of the negatively charged hyperons is substantially suppressed with re- 
spect to that in neutron-star matter, since charge neutrality is handled by the electrons, 
which must maintain a sizable fraction throughout the entire density range. 

2. The formation of the S~ gets delayed, and thus the A is the most favorable hyperon to 
be produced at the lower densities. The lack of competition from the S~ helps the A, 
precisely as for the neutron star matter case when the S's were excluded by assuming 
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repulsive S-nm interactions. As a result, we find that A's appear first in the matter, 
at slightly lower densities than when produced in neutron star matter, and it is the 
S~ which comes in a close second. A similar effect governs the fractions of the E*^, T,~^ 
and H*^ with respect to the S~. 

3. The total hyperon fraction continues to be smaller than in neutron star matter through 
the enitre density range, although the A's still accumulate to a considerable fraction. 

4. The effect of hyperon formation on the lepton population is also worth attention. Once 
again, the formation of hyperons makes possible the reduction of the neutron fraction 
without a coupled growth in the proton fraction. Thus the difference /i„ — /ip decreases, 
forcing a similar trend in /ig — fiue- The immediate result is that hyperon formation 
is closely followed by a reverse of the trend in the electron and neutrino populations, 
with respect to those of nuclear matter. The neutrino population begins to rise, and 
quickly overtakes its value at nuclear saturation density of about 0.06 per baryon, and 
even crosses 0.1 per baryon at high densities. 

Again, the results found here with the effective EoS bear high resemblance to those reported 
in works which used more sophisticated models 0. 

The above results depend only little on temperature, and similar equilibrium compositions 
are found for any significant temperature (T > 10 MeV). The temperature determines the 
precise hyperon populations at low densities, but these are insignificantly low anyway. There 
is some model dependence with respect to the parameterization of the baryon-baryon effective 
interaction (the values of 7 and 6), in a similar fashion as for neutron star matter. Again, 
very little dependence is found regarding the value of the hyperon-hyperon well depth, which 
is of even less importance in the present the matter is even more nucleon dominated. 

Similar results are found for the case of a repulsive S-nm interaction, excluding S's from 
the matter. The equilibrium compositions for proto-neutron-star matter with 7 = 5 = | 
and —Vh^h2 = 40 MeV, but when the S's are excluded, are plotted in Fig. 10c. Since the 
formation of the S~ is suppressed even when allowed, A's dominate inherently and the A 
fraction is rather insensitive to the presence of S's. Excluding them does create an increase, 
obviously, in the S fraction. 

The outstanding point of proto-neutron star composition is that it generates a different 
equation of state than for neutron star (neutrino-free) matter. In nuclear matter the effect 
is trivial: the greater symmetry among the nucleons lowers the total Fermi energies and 
the symmetry potential energy density. Hence, equilibrium nuclear matter with trapped 
neutrinos is inherently softer than neutrino-free nuclear matter for the same nuclear-force 
model. 



24 



The picture is quite the opposite when hyperons enter the game. As discussed above, 
the more symmetric nuclear fraction somewhat suppresses hyperon formation in neutrino- 
trapped matter. Since hyperon formation imposes a softer equation of state by lowering 
kinetic and potential energies, it is clear that when hyperons are allowed to appear the equa- 
tion of state for neutrino-trapped matter, which has a smaller hyperon fraction will be stiffer 
than the EoS for neutrino free matter. This principle was pointed out in [7,8], and is laid out 
in detail in [|]. The equations of state calculated here do indeed reproduce these expected 
results. Fig. 11 compares pressure vs. baryon density (which is the relevant quantity when 
using the same model for different thermodynamic conditions), for 5 = 7 = | and —Vhj^h2 =40 
MeV. As is expected for nuclear matter (Fig. 11a) the neutrino-trapped case (dashed line) 
is clearly softer (except at very low densities, where the finite-temperature lepton pressure is 
of significance). Matter with hyperons presents the opposite state of affairs (Figs, lib and 
11c): as is clearly seen, the neutrino-trapped cases do indeed have stiffer equations of state 
than the neutrino- free ones (solid lines). 



6 Results: Respective Neutron Star Masses 

A natural test for the constructed equations of state is the calculation of the neutron star 
properties they predict. Most important of these is, of course, the mass/ cent eral- density 
relation: the mass sequence M{ec), of the neutron star for a given equation of state, and 
especially the maximum mass predicted by the equation. This maximum mass must be 
at least 1.4 Mq, which is the well determined observed mass of the Hulse- Taylor pulsar 
(see 1^ for a recent review of observed and calculated neutron star properties), where 
= 1.989 X 10^^ gm is the standard notation for one solar mass. 

In this study we limit our discussion to calculating the static neutron star masses, noting 



that rotation may add as much as 0.1-0.2 Mq [^. The static masses are calculated by 
numerical integration of the Tolman-Oppenheimer-Volkoff equations [^] of hydrostatical 
equilibrium, namely (setting G=c=l, so that all physical quantities are given in centimeters, 
Mq = 1.477 X 10^ cm): 

^ = 4-^- ; li7 = -7r(i + T)(i + — 1(1--) (32) 

(and — : 
dr 
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These equations are solved following the recipe of Arnett and Bowers by setting a central 



density £c in the star and integrating outwards until reaching zero pressure. The equation of 
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state of Baym, Pethick and Sutherland |^ is used for the subnuclear densities. The specific 
choice of the EoS for these densities is actually unimportant for the mass sequence, as the 
fraction of the neutron star mass originating from the subnuclear densities is insignificant. 



6.1 Neutron Star Masses 

The predicted mass sequences for the three equations of state for nuclear matter (5 = 2, | 
and |) are plotted in Fig. 12 as a function of the central density. As is expected from 
the large difference between these equations (Fig. 6), the neutron star mass sequences vary 
considerably from one equation to another. For the stiffest equation (5 = 2), at ec = l-25x 10^^ 
gm/cm^, after which the EoS breaks causality, the mass begins to saturate at the very high 
value of 2.6 Mq. The softest equation, which is still stiff in comparison to most of the more 
sophisticated nuclear matter models, yields a maximum mass of about 2.2 Mq. 

The consequences of the relative similarities between the equations of state for matter with 
hyperons are even enhanced in the neutron star mass calculations. As noted in Section 5, the 
dominant factor for the EoS is the number of species that are available for production in the 
matter, the case of a repulsive S-nm interaction being clearly different than the attractive 
one, and both are very different than the nuclear matter EoS. 

Fig. 13 presents six plots of neutron star mass sequences. In all cases the hyperon-hyperon 
well depth is chosen as —Vh^h2 = 40 MeV, and the curves differ in the short range repulsive 
forces model, i.e. in the values of 7 and 6. For comparison, the mass function of the softest 
nucler matter equation 5 = | is also plotted. The outstanding result is that while in each case 
the masses corresponding to low densities are very model dependent in terms of the short 
range interactions (the values of 7 and 6), the maximum masses are very similar for different 
equations of the same type (including or excluding the S's). The behavior at low central 
densities must be different, as it results mainly from the nucleon-nucleon interactions, which 
determine the precise model dependence. The similarity of the mass sequences at high central 
densities is a direct consequence of the hyperon "pressure control" mechanism, discussed in 
the previous Section. This mechanism is basically dependent only on the availability of the 
various hyperon species for formation in the matter - and hence, the clear distinction between 
matter which may contain S's and matter which does not. 

These results clearly indicate that the appearance of hyperons makes the maximum mass of 
neutron stars sensitive only to the gross properties of the baryon-baryon interactions, and 
almost insensitive to specific details. While these results generated with effective equations 
of state can hardly be considered as proof of such a trend, we do believe that they do serve 
as indication that it is the number of available species that dominates the neutron star mass 
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function at high density. 

The precise values of the resulting maximum masses are also of interest. Since the baryonic 
matter EoS is softer than that for nuclear matter, the maximum neutron star masses are 
expectedly smaller. The maximum mass in the various models for matter without S's are 
about 1.6-1.7 Mq, while the softer matter with E's holds maximal masses of only about 
1.4-1.5 Mq. The proximity of the latter value to the upper limit of observed neutron star 
masses (1.44 M0 of PSR 1913+16) may be viewed in two ways: as an indication that S's are 
bound in the matter, thus limiting the maximal mass to less than 1.5 Mq; or on the contrary 
- since the EoS with S's yields such a border-line result, then perhaps it is more plausible 
that E's are unbound in the matter. This possibilty is enhanced by the recent suggestions 
that the mass of the Vela Pulsar is above 1.55 Mq, based on analysis of its X-ray source 



properties [EC 



While we are unaware of an explicit indication to a hyperon induced "pressure control" mech- 
anism in previous works, it is interesting to note that most studies which include hyperons 
in neutron star matter find that the maximum masses lie within a narrow range, typically 
1.4— 1.7Mq (see Refs. [|], ||, |, |10|, |3^ and other references therein), not unhke the results 



presented here. On the other hand, published neutron star maximum masses for nuclear 
equations of state do indeed vary over a larger range (1.5 — 2.5 Mq), and are usually larger 
than 1.8 Mq [ibid.]. Both these results are consistent with the hyperon induced "pressure 
control" mechanism discussed above. 

We complete this study by remarking that all results show little dependence on the hyperon- 
hyperon well depth, as is expected for nucleon-dominated matter. Choosing a large well 
depth along with a steep high density repulsion (7 = 2) creates greater repulsion at high 
densities, leading to a stiffer EoS and a slightly higher maximal mass. However, we find 
that this effect changes the maximum mass by AM < O.IMq between the case of strong 
hyperon-hyperon interactions and the case when these interactions are turned off. 



6.2 Proto Neutron Star Masses 

Realisitic proto-neutron star calculations require consistent thermodynamic solutions, in 
order to establish the physical temperature profile for the star (as the star's physical state 
is probably closer to constant entropy per baryon rather than constant temperature). We 
thus limit ourselves to simplified static solutions, in order to recover the key feature of 
proto-neutron star masses. 

Since the EoS for nuclear matter stiffens when neutrinos diffuse out of the proto-neutron 
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star, the maximal mass of the fully evolved neutron star will be greater than that of the 
proto-neutron star. Hence, if a stable proto-neutron star is formed in the supernova (and 
does not collapse initially into a black hole), it will also have a stable configuration in the 
fully evovled state. 

The opposite trend in matter with hyperons also brings about the possibility of the reverse 
physical scenario: since the equation of state for matter with hyperons softens as the neu- 
trinos diffuse outward (and the hyperon population increases), the maximal mass of the 
proto-neutron star will be greater than that of the fully evovled star, under the same physi- 
cal assumptions. Thus, the stage is set for a delayed coUpase of the proto-neutron star into 
a black hole, if its mass will be too large for the neutrino-free equation of state. Such a 
scenario has been motivated by the combination of a neutrino pulse in SN1987A (so there 
was no initial collapse into a black hole), and the fact that no pulsar has been detected; see 
Refs. [7-10]. 

We reproduced the necessary features described above in the effective EoS scheme, when 
deriving mass sequences for proto-neutron stars and evolved neutron stars (Fig. 11). For 
nuclear matter the EoS for proto-neutron star (neutrino-trapped) is softer than for evolved 
neutron star (neutrino-free) and the maximal mass drops from ~ 2.4 Mq to ~ 2.0 Mq for 
the proto-neutron star. Since the EoS for neutrino-trapped matter is stiffer for matter with 
hyperons, the maximal mass of the proto-neutron star is found to be larger for the hyperonic 
case. For matter with all types of hyperons the maximum mass of the proto-neutron star is 
almost 1.7 Mq (with respect to a httle over IAMq for the evolved neutron star). Matter with 
no E's is less affected by the neutrino diffusion (as explained above, the E's are somewhat 
suppressed in neutrino-trapped matter) but the proto-neutron star still has a maximum mass 
of about O.IM0 higher than that of the evovled star (l.SM© vs. 1.7 Mq). This trend appears 
for the other various cases of assumptions made about the baryon-baryon interaction and 
shows only little model dependence, as in the case of the evolved neutron star masses. 

7 Conclusions and Discussion 

The effective equation of state described in the present study is designed to allow rapid cal- 
culations and extensive parameter surveys of high density matter rich with strange baryons. 
This equation generalizes the Lattimer-Swesty EoS for nuclear matter, and is constructed so 
that it will reproduce the physical features attributed to such matter, while not presuming 
to correctly describe the microscopic details of the underlying physics. Its construction is 
performed by basing the EoS on local (effective) potentials for the various componenets of 
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the baryon-baryon interaction, which are constrained as much as possible by experimental 
data from hypernuclei. 

In view of the large uncertainties in the baryon-baryon interactions, especially at high densi- 
ties, we examine the properties of the EoS with various assumptions regarding these interac- 
tions. We investigate sensitivity to the density dependence of the repulsive term and remark 
on the unknown power of the hyperon-hyperon interactions. Special attention is given to the 
consequences of a recent analysis of E-atoms data, which suggests a strong isoscalar repul- 
sion in the E-nuclear matter (E-nm) interaction (very different than the commonly assumed 
interaction), and excludes E hyperons from forming in the high density matter. 

We find that the effective equation of state yields similar equilibrium compositions and 
equations of state as found in equivalent RMF models. The main features of the hyperon 
formation in zero-temperature neutrino free matter (typical of evolved neutron star interiors) 
are: 

1. Hyperons begin to appear in the matter at about twice the nuclear saturation density 
Ps ^0.155 fm~^, and at 2.5ps the hyperon fraction in the matter is already significant. 
First to appear are the (due to its negative charge) and the A, or the A alone, 
depending on the assumed E-nm interaction. Other hyperons follow at higher densities. 

2. Hyperon production is followed by considerable deleptonization, as negatively charged 
hyperons replace electrons and muons in maintaining charge neutrality. 

We find that these features are typical of all models, although the details arc model- 
dependent as far as the density dependence of the baryon-baryon repulsive interaction terms. 
Since the matter is dominated by nuclcons in most of the density range in question, there is 
little dependence on the hyperon-hyperon interactions. 

The key effect of the appearance of hyperons in the matter is the softening of the equation of 
state, due to the lowering of both repulsion and Fermi energies as the total baryon density is 
divided among a larger number of species. We emphasize that hyperons serve as a "pressure 
control" mechanism in the high density matter, since the extent of their appearance is directly 
related to the details of the nuclear interactions. Thus, while for nuclear matter a more 
powerful nucleon-nucleon repulsion creates a stiffer EoS, in baryonic matter an enhanced 
production of hyperons will occur, restoring the gross properties of the EoS. The immediate 
result of this "pressure control" by the hyperons is that all EoS including hyperons arc very 
similar, the only significant model dependence being on the number of available species (i.e. 
the general features of the E-nm interaction discussed above) . 
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The most notable manifestation of this effect is in the predicted maximum masses of neutron 
stars. By solving the static mass functions of neutron stars for the various equations, we 
find that the hyperon-imposed similarity in the equations of state leads to very moderate 
specific-model-dependence of the maximum masses. While the maximum mass for nuclear 
matter is highly model dependent, the maximum masses for matter with hyperons depend 
almost entirely just on the number of available species. Specifically, wc find that for matter 
which includes S's, the maximum mass for all equations of state are about 1.4— l.SM©, while 
for the case of excluded S's, maximum masses are all in the range 1.6— 1.7Mq. Such analysis 
with an effective equation of state may hardly be considered as a proof of this point, but we 
do believe that the general trend described above indicates a genuine physical principle. 

It is interesting that the result for matter with E's yields just about the maximum observed 
neutron star mass (the small difference can be bridged by allowing also for rotation and 
finite temperature). This may be viewed in two ways: one could either argue that this is too 
border-line a result, implying E's should be excluded, or rather find this an indication that 
the maximum mass is indeed close to the observed value, as determined by the hyperons. 
Clearly further analysis is needed to elaborate on this point. 

We briefly investigate the properties of proto-neutron star matter, which has a finite tem- 
perature and in which neutrinos are assumed to be trapped, so that the total lepton number 
remains constant. Such matter is relevant to the simulation of proto-neutron star evolu- 
tion and core-collapse process in supernovae. Again we find that hyperons appear at about 
2ps and accumulate significantly by 2.5ps, with very weak model dependence. Hyperon ap- 
pearance changes the trend of the neutrino population, which begins to increase rather than 
decrease at higher densities, since the hyperons indirectly lower the electron fraction required 
to maintain charge neutrality. These results are in good agreement with studies performed 
with more sophisticated models. 

The formation of hyperons softens the equation of state for core-collapse matter as well, and 
once again serves as "pressure control" for the nucleons. The softening occurs, of course, 
at the density of hyperon accumulation, which is border-line with regard to the densities 
reached during the supernova. None the less, a softer equation of state is appealing, as it 
will strengthen the shock and the neutrino delayed shock, thus strenghtening the explosion. 
In the proto-neutron star case such a softening is clearly important, and we do reproduce the 
well-known tendency of the maximum mass for matter with trapped neutrinos being larger 
than that for neutrino free matter. This feature, which is typical for matter with hyperons 
and opposite to the one prevalent in nuclear matter, opens the way for the "delayed collapse" 
scenario, motivated by the lack of an observed pulsar in the Supernova 1987A remnant. In 
such a scenario, a newly formed neutron star could exist in a meta-stable configuration. 
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collapsing into a black hole when the neutrinos diffuse out of the star leaving it with too 
large a mass. 

All of the above results serve as indication that the effective equation of state may be treated 
as a trustworthy tool in the analysis of high density matter with strangeness. Nonetheless, 
some further extension of the equation should be considered. One area for improvement is 
the description of the high density behavior of the forces. A possible approach would be to 
add another repulsive term c'p'^', where 7' will be greater than 7 (which could then be set 
to a low value allowing consistently for a low nuclear incompressibility) . This new repulsive 
term could be set universaly for all species (e.g, if it may be attributed to quark degrees of 
freedom). This is somewhat similar to a suggestion by Lattimer, Swesty and Myra ||2^, that 



the repulsive term in the energy density formula be amended to cp^'^^ j (1 + dp^ as an ad 
hoc method to artificially soften the equation. 

Other possible improvements are introducing finite ranges to the forces, and a more com- 
prehensive treatment of effective masses. Clearly, any future experimental and theoretical 
constraints on the nuclear, nucleon-hyperon and hyperon-hyperon forces must be incorpo- 
rated into the equation and its parameters as well. 

Introduction of vr and K condensates should also be considered, although viewed as unlikely 
by several current works. The same applies to the formation of A isobars, which are found 
to be excluded from the matter by using standard assumptions regarding its interactions and 
mass. A natural extension of the equation would be to include a deconfined quark phase at 
high densities, for which an effective equation of state is quite simple p3 . 



Further work is also required regarding the implications of strangeness on high density mat- 
ter in astrophysical contexts. Up to date no neutron star property that can be uniquely 
connected to hyperon formation has been suggested, although observed cooling rates have 
been shown to place some limits on the matter composition (a large hyperon fraction along 
with a significant lepton population would produce much too large cooling rates through the 
direct Urea mechanism). Studies of possible implications of hyperon formation on rotational 
properties of neutron stars, along with magnetic, superconducting and superfiuid behavior of 
such matter, and crust properties with relation to observed glitches should also be considered 
(far more extensive studies on these subjects have been carried out for strange-quark-stars, 
following the still debated stable-strange- matter hypothesis). We note that a recent work by 



Glendenning, Pei and Weber |^ suggests that the existence of a quark-matter core (which 
will definitely include a non- zero-strangeness-fraction) could be identified by its effect on the 
pulsar breaking index, if such a core builds up during the pulsar spin down. 

A unique signature of hyperon (and generally strange particles) production in proto-neutron 
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star would also be of considerable interest, especially regarding the "delayed-collapse" sce- 
nario and whether it may lead to an appreciable effect on the neutrino pulse from the 
proto-neutron star [Q. Studies along these lines have begun in recent years, and are well 
worth continuing. 

Last but clearly not least, the results presented in this study imply that hyperon formation 
may be of importance in core-collapse supernovae, if densities equal or greater to twice 
the nuclear saturation density are reached in the collapsing core prior to the "bounce". A 
softened equation of state might claim its role in the shock- wave/neutrino-delayed-shock 
process, presumably strengthening the explosion. 
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baryon-baryon interactions. We further wish to thank Eliahu Friedman, Itamar Lichtenstadt, 
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Figure Captions 



Fig. 1. The density dependence of A potential in nuclear matter, Van{pn) = o.knPn-^cxnPni 
for different values of the exponent 7 in the repulsion term, using the parametrs sets of MDG 
[ p^ given in Table 2. The arrow marks the nuclear saturation density, ps = 0.155 fm~^. 

Fig. 2. The equilibrium compositions of nuclear matter (n-p-e-yu) as a function of baryon 
density for the case 5 = 2; see Eq. (j^) and Table 1. Other models of the nuclear interaction 
(different values of 5) yield similar equilibrium compositions. 

Fig. 3a. The equilibrium compositions for matter containing hyperons as well as nucleons 
and leptons, for 5 = 7 = 2 (see Eq. (|D), and for the hyperon-hyperon well depth —Vh^h2 = 40 
MeV. 

Fig. 3b. Same as 3a but with 5 = 7= |. 
Fig. 3c. Same as 3a but with 5 = 7= |. 

Fig. 4a. The equilibrium compositions for matter containing A's and H's as well as nucleons 
and leptons, but no S's (when a repulsive E-nm interaction is assumed), for 5 = 7 = 2, and 
a hyperon-hyperon well depth —Vhj^h2 = 40 MeV. 

Fig. 4b. Same as 4a but with 5 = 7= |. 

Fig. 4c. Same as 4a but with 5 = 7= |. 

Fig. 5. Same as 3a but with hyperon-hyperon interactions "turned off" . 

Fig. 6. Equations of state for nuclear matter (n-p-e-/i) for the different models of the nucleon 
potential energy density. The equations for 6 = 2 and | break causality at 1.25 x 10^^ gm/cm^ 
and 1.94 x 10^^ gm/cm^, respectively, as indicated by the dashed curves. 

Fig. 7. Equations of state for high density matter with nucleons and leptons alone; nucleons, 
hyperons and leptons; and A's and S's as well as nucleons and leptons, but no S's (for a 
repulsive E-nm interaction). The parameters used are 5 = 7= |, and —Vhih2 = 40 MeV. 

Fig. 8. Equations of state for high density matter with nucleons, hyperons (assuming an 
attractive E-nm interaction) and leptons, for different choices of the baryon potential energy 
density dependence, 5 = 7 = 2, | and |, and a hyperon-hyperon well depth — V/i^/ij = 40 MeV. 
The equations of state correspond to the equilibrium compositions of Figs. 3a, 3b and 3c. 

Fig. 9. The effect of hyperon formation on the adiabatic index F of the EoS for the case of 
5 = 7= |, and —Vh-^h2 = 40 MeV, compared to the corresponding EoS for nuclear matter. 
Each "dip" in the curve of the adiabatic index corresponds to the appearence of a new 
hyperon species. 
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Fig. 10a. The equilibrium compositions of neutrino trapped nuclear matter (n-p-e-i/g) as 
a function of baryon density for the case S — ^, with a lepton fraction oi Yi — 0.36 and 
temperature T — 20 MeV. The equihbrium compositions for other models of the nuclear 
interaction (values of S) are very similar. 

Fig. 10b. The equilibrium compositions for neutrino-trapped matter containing hyperons as 
well as nucleons and leptons, where 5 = 7= |, and —Vhih2 = 40 MeV. 

Fig. 10c. The equilibrium compositions for neutrino-trapped matter containing hyperons as 
well as nucleons and leptons, but no E's (for a repulsive E-nm interaction), for 5 = 7 = | and 
-Vk,h, = 40 MeV. 

Fig. 11. Equations of state for neutrino- free matter (solid lines) and neutrino-trapped matter 
(dahed lines) for 5 = 7 = |, for the different versions of high density matter: nucleons only 
(a), nucleons and all hyperons (b) and nucleons, A's and S's but no E's (c). It is noticeable 
that the EoS for the neutrino-trapped case is softer for nuclear matter, but stiffer for matter 
with hyperons. 

Fig. 12. Neutron star mass vs. central density mass-sequences for nuclear matter (n-p-e-//) 
for the different models of the nucleon potential energy density. 

Fig. 13. Neutron stars mass vs. central density mass-sequences for high density matter with 
all hyperon species and with only A's and S's but no E's (for a repulsive E-nm interaction), 
for different values of 7 and 6. In all cases the hyperon-hyperon well depth is —Vh^h2 = 40 
MeV. Also plotted is the neutron star mass sequence for the softest equation for nuclear 
matter 5 = |. 
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Table 1: 

Values for the parameters of the terms corresponding to the nuclear potential energy density 
of Eq. (H), as derived from the Lattimer-Swesty equation of state with m* = m. 



6 


ttNN 


bNN 


cnn 




[MeV fm^] 


[MeV fm^] 


[MeV fm^^] 


2 


-784.4 


214.2 


1936.0 


5 
3 


-935.4 


214.2 


1557.2 


4 

3 


-1384.6 


214.2 


1672.8 



Table 2: 

Values for the parameters of the terms corresponding to the A-nm interactions in the local 
potential energy model of Eq. for a well-depth of —Van = '^^ MeV at nuclear saturation 
density. 



7 




Can 




[MeV fm^] 


[MeV fm^^] 


2 


-340.0 


1087.5 


5 
3 


-387.0 


738.8 


4 
3 


-505.2 


605.5 
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Table 3: 

Values for the parameters of the terms corresponding to the S-nm interactions in the local 
potential energy model of Eq. (1), for a well-depth of — VsAr = 28 MeV at nuclear saturation 
density, assuming a E-nm interaction equal to the vVnm one, except for a nonzero isovector 
term twice that of the N-nm interaction. 



7 










[MeV fm^] 


[MeV fm^] 


[MeV fm^T] 


2 


-340.0 


214.2 


1087.5 


5 
3 


-387.0 


214.2 


738.8 


4 

3 


-505.2 


214.2 


605.5 



Table 4: 

Values for the parameters of the terms corresponding to the S-nm interactions in the local 
potential energy model of Eq. (|^), accepting the strong E-nm isoscalar repulsion argued 
for in Ref. |]TB[. The b and B parameters are given in fm. 



a 


bo 


Bo 


bi B, 


0.5 


2.1 


-4.3 


-1.0 



Table 5: 

Values for the parameters of the terms corresponding to the S-nm interactions in the local 
potential energy model of Eq. (^, for a well-depth of — VHAr = 24 MeV at nuclear saturation 
density. 



7 




bsN 






[MeV fm^] 


[MeV fm^] 


[MeV fm^T] 


2 


-291.5 





932.5 


5 
3 


-331.8 





663.4 


4 

3 


-434.4 





520.1 
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Table 6: 

Values for the parameters of the terms corresponding to the various hyperon-hyperon in- 
teractions in the local potential energy model of Eq. assuming the hyperon-hyperon 
well-depth at a density equal to the nuclear saturatrion density to be — V/^/ij = 40 MeV. 
Only the E's are assumed to have a nonzero isovector component in their interactions. 



7 




bAY=bEY 




cyy 




[MeV fm^] 


[MeV fm^] 


[MeV fm^] 


[MeV fm^T] 


2 


-486.2 





428.4 


1553.6 


5 
3 


-552.6 





428.4 


1055.4 


4 
3 


-723.2 





428.4 


869.0 
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figure 5 
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